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I.  SUMMARY 


The  combined  system  consisting  of  the  baseline  flexible  structure  modified  by  the 
system  of  active  controllers  is  considered  as  a  unified  dynamical  system.  Techniques 
based  on  computer  algebra  (MACSYMA)  are  used  to  derive  expressions  for  the 
transfer  functions  of  the  modified  system,  using  the  known  transfer  functions  of  the 
baseline  flexible  structure  and  the  feedbacK  gains  of  the  active  controller.  The  roots  of 
the  characteristic  polynomial  of  this  transfer  function  give  the  system  resonant 
frequencies  and  damping  parameters.  Using  the  computer  algebraic  system 
MACSYMA,  expressions  for  these  parameters  which  are  explicitly  dependent  on  the 
output  feedback  gains  of  the  active  controller,  are  presented.  For  lightly  coupled 
modes,  simple  relations  are  obtained  between  the  modal  parameters  and  the 
coordinates  of  the  sensor/actuator  pairs  as  well  as  the  displacement  and  velocity 
feedback  gains.  These  results  permit  the  parametric  study  of  the  placement  of  the 
resonant  frequencies  and  damping  parameters  of  the  combined  system,  as  functions 
of  the  feedback  gains.  Numerical  examples  are  used  to  illustrate  the  application  of 
these  results  to  the  calculation  of  active  controller  feedback  gains  based  on  the 
requirement  that  certain  modes  have  specified  modal  damping  while  the  closed-ioop 
frequencies  remain  unchanged.  Also,  using  the  example  of  a  cantilevered  uniform 
beam  controlled  by  a  single  sensor/actuator  pair,  numerical  results  are  used  to 
illustrate  the  sensitivity  of  the  closed-loop  modal  parameters  to  the  placement  of  the 
sensor/actuator  pair  as  well  as  the  feedback  gains.  Such  results  help  to  answer 
questions  about  optimal  placement  of  sensor/actuator  pairs  for  the  active  control  of  a 
flexible  structure.  [Key  words:  Active  Control;  Smart  Structures;  Computer 
Algebra;  Dynamics  of  Flexible  Structures] 
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II.  NOMENCLATURE 


aj,  aj 

= 

measurement  coordinates  for  ith  and  jth  controllers 

bj.bj 

= 

force  application  coordinates  for  ith  and  jlh  controllers 

CjO) 

= 

symbolic  coefficient  of  j*h  power  of  p  in  the  i-term  derivation 

G(x,4,t,x) 

= 

Green’s  function 

9j-  9j 

= 

displacement  feedback  gains  of  ith  and  j,h  controllers 

hi.hj 

= 

velocity  feedback  gains  of  i^  and  jth  controllers 

I 

= 

identity  matrix 

L 

= 

number  of  discrete  attachments  (controllers) 

*-x>t{ ) 

= 

partial  differential  operator 

P 

= 

Laplace  variable 

POk 

s 

kth  parameter  of  baseline  system 

Q(x,t),  Q(x,p) 

= 

system  response  and  its  Laplace  transform 

t 

= 

time 

W(x£,p) 

= 

system  transfer  function 

Wci(p) 

= 

transfer  function  of  ith  controller 

w(x,t),  w(x.p) 

= 

forcing  function  and  its  Laplace  transform 

x,  X1 ,  X2 

= 

spatial  coordinates 

= 

coefficient  of  the  i*h  power  of  p  in  characteristic  polynomial 

P>  Pk 

= 

modal  weights 

2 


50 

{tf.Yj 

n 

^Ok 


a>a0ni’°1  ni 

t 


w.WQni.^i  ni 


Dirac  delta  function 

vector,  with  elements  defined  in  Equation(15) 
spatial  coordinate 

k*h  orthonormal  modal  function  for  baseline  system 

exponential  growth  rate 

time 

frequency 


{to},  COj 


[O],  Qjj 


vector,  with  elements  defined  in  Equation(13) 
spatial  coordinate 

matrix,  with  elements  defined  in  Equation(12) 


III.  INTRODUCTION 

This  research  effort  has  been  concerned  with  the  development  of  new  techniques  for 
the  analysis  of  the  dynamic  response  of  distributed  parameter  flexible  structures  which 
are  being  controlled  actively  at  discrete  locations.  This  effort  is  a  collaboration 
between  Clark  Atlanta  University  and  AEDAR  Corporation.  This  work  grew  out  of  the 
advancements  made  by  AEDAR  in  the  application  of  computer-algebraic  techniques  to 
the  Green's  function  analysis  of  structural  systems  modified  by  discrete  substructures 
attached  to  them  (Fabunmi,  1989). 

The  technology  of  active  control  of  structures  has  been  the  subject  of  considerable 
research  interest  over  the  past  decade  in  part  because  of  the  need  to  suppress 
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excessive  vibrations  associated  with  the  deployment  of  large  flexible  structures  in 
space  (Various  Authors,  1986;  Atluri  and  Amos  (Ed.),  1988).  Because  space-borne 
structures  cannot  afford  tne  weight  penalties  of  classical  vibration  control  devices  such 
as  absorbers  or  isolators,  a  lot  of  effort  has  been  devoted  to  various  means  of  actively 
controlling  the  dynamic  characteristics  of  these  structures.  These  techniques  use  an 
external  source  of  energy  to  apply  controlling  forces  (and/or  moments)  on  the  structure 
which  are  determined  in  some  relationship  to  the  measured  or  estimated  response  of 
the  structure.  More  advanced  implementations  of  active  structural  control  involve  the 
embedding  of  sensors,  actuators  and  processors  in  the  structure  itself.  Such  "smart 
structures"  are  able  to  adjust  the  characteristics  of  their  controllers  e  g.  feedback  gains, 
in  response  to  changing  dynamical  environments. 

The  most  popular  approach  to  the  design  of  the  active  control  schemes  follows  the 
paths  of  modern  control  theory  which  involves  optimal  state-space  feedback  control 
(Various  Authors,  1984,  1986;  Atluri  and  Amos  (Ed.),  1988;  O’Donoghue  and  Aiiuri, 
1985;  Horner  and  Walz,  1985),  or  output  feedback  control  (Meirovitch,  1988;  Garcia 
and  Inman,  1990).  A  finite  order  mathematical  model  of  the  structure  is  required.  The 
state  variables  are  the  [generalized]  displacements  and  [generalized]  velocities  of  the 
structure.  In  the  case  of  state-space  feedback,  the  state  of  the  system  is  estimated  from 
measurements  at  selected  coordinates  and  this  estimate  is  used  to  derive  the 
feedback  gains,  using  a  method  based  on  Pointryagin's  principle  for  solving  a 
constrained  optimization  problem.  This  method  involves  the  computation  of  a  positive- 
definite  matrix  satisfying  the  algebraic  matrix-Riccatti  equation  (Junkins  and  Rew, 
1988).  Output  feedback  control  does  not  use  the  entire  state-space  estimate  for 
feedback;  instead  only  the  measured  responses  are  used.  The  advantage  is  that  the 
practical  implementation  of  the  controller  is  simpler  and  errors  associated  with  the 
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estimation  of  unmeasured  responses  are  eliminated  (Garcia  and  Inman,  1990).  For 
designs  based  on  these  methods,  the  control  strategy  is  specified  in  terms  of  the 
minimization  of  an  objective  function.  For  a  given  objective  function  and  a  set  of  initial 
conditions,  the  controller  feedback  gains  are  calculated  once  and  for  all.  A  question 
that  is  often  asked  is  whether  or  not  optimal  control  necessarily  implies  intelligent 
control.  For  the  control  to  be  considered  intelligent,  the  system  must  have  the  ability  to 
alter  its  control  strategies,  and  readjust  the  feedback  gains,  based  on  some  reasoning. 
An  example  of  such  reasoning  is  to  have  a  controlled  structure  which,  upon  sensing 
(using  its  embedded  sensors)  that  the  spectral  distribution  of  its  current  external 
excitation  is  close  to  0^9  or  more  of  its  resonant  modes  (the  parameters  of  which  has 
been  stored  in  the  memory  of  its  embedded  processor),  can  alter  the  strategy  of  its 
active  controller  such  that  the  necessary  feedback  gains  are  computed  (using  its 
embedded  processor)  which  will  maximize  the  modal  damping  of  the  most  highly 
excited  modes.  This  type  of  approach  requires  that  algorithms  be  available  for  directly 
calculating  the  controller  feedback  gains,  based  on  specifications  of  required  closed- 
loop  modal  parameters.  The  development  of  such  algorithms  has  been  one  of  the 
motivations  of  this  research  effort. 

Recent  advances  in  computer  algebra  have  made  available  symbolic  manipulation 
facilities  which  extend  the  tools  of  algebra  and  integro-differential  calculus  beyond  the 
traditional  limits  (Rand,  1984;  Pavelle,  1985).  Techniques  based  on  computer  algebra 
were  reported  by  Fabunmi  (1989),  which  permit  the  derivation  of  the  transfer  functions 
(Laplace  transform  of  the  Green's  Functions)  of  the  system  resulting  from  the 
attachment  of  discrete  dynamic  substructures  to  a  distributed  parameter  base-line 
structure.  It  is  assumed  that  the  algebraic  forms  of  the  transfer  function  of  the  base-line 
structure  as  well  as  those  of  the  discrete  attachments  are  known.  The  mathematical 
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form  chosen  for  system  transfer  functions  permits  the  direct  determination  of  the 
system  parameters  as  the  complex  values  of  the  Laplace  variable  at  whicf 
singularities  of  the  transfer  function  occur.  For  the  class  of  controllers  where  the 
measured  output  are  the  feedback  variables  such  as  displacements  and  velocities, 
the  resulting  system  is  mathematically  equivalent  to  that  of  the  attachment  of  discrete 
"substructures”,  the  transfer  functions  of  which  are  given  by  expressions  involving  the 
gain  constants  and  the  Laplace  variable. 

This  report  has  been  organized  as  follows.  In  the  section  IV,  the  equations  which  were 
used  to  derive  the  effects  of  active  controllers  on  the  system  parameters  are  derived. 
This  section  includes  some  pertinent  material  from  (Fabunmi,  1989),  for 
completeness.  Section  V  presents  the  expressions  for  the  characteristic  polynomials 
of  the  closed-loop  system,  obtained  using  MACSYMA  on  the  Symbolics  3620 
workstation.  This  is  followed  by  examples  of  how  these  results  can  be  used  to 
calculate  the  controller  feedback  gains  based  on  specified  parameters  of  the  closed- 
loop  system  (section  VI),  and  in  section  VII,  using  the  example  of  a  cantilevered 
urniarm  uearn  controlled  by  a  single  Se.iSunactuator  paii,  numerical  results  are  used  to 
illustrate  the  sensitivity  of  the  closed-loop  modal  parameters  to  the  placement  of  the 
sensor/actuator  pair  as  well  as  the  feedback  gains.  The  basic  conclusions  of  this  report 
which  are  presented  in  the  section  VIII,  are  that  (1)  new  tools  based  on  computer 
algebra  have  been  developed,  for  the  analysis  of  system  characteristics  of  actively 
controlled  structures,  (2)  alternative  techniques  for  the  design  of  active  controllers 
have  been  presented,  which  make  it  possible  to  design  an  adaptive  controller  for 
which  different  schedules  of  feedback  gains  can  be  used  to  adjust  the  system 
parameters  as  needed,  in  order  to  minimize  dynamic  response  to  external  excitations, 
and  (3)  these  results  make  it  possible  to  select  locations  for  sensors  and  actuators,  at 
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which  the  closed-loop  system  parameters  have  the  desired  sensitivities  to  the  velocity 
and  displacement  feedback  gains. 

IV.  ANALYSIS  OF  ACTIVE  CONTROLLER  EFFECTS 

The  objective  of  this  section  is  to  present  the  formulation  of  the  equations  that  were 
used  to  derive  the  transfer  function  of  the  system  resulting  from  the  attachment  of  a 
finite  number  of  discrete  linear  output  feedback  controllers  to  a  distributed  parameter, 
baseline  system  such  as  a  flexible  structure.  These  derivations  follow  the  same  lines 
as  those  presented  in  (Fabunmi,  1989).  The  dynamic  response  of  a  distributed 
parameter  system  are  solutions  to  partial  integro-differential  equations  which  can  be 
represented  operationally  as: 

^„{Q(*  2,t)}  =  w(x11t)  (1) 

where  Lxt  {  }  is  an  integro-differential  operator  which  maps  the  responses  Q(x2,t)  on 
to  the  excitations  w(x1,t)  subject  to  appropriate  boundary  and  initial  conditions  on 
Q(x2,t);  x1  is  in  the  spatial  domain  of  the  excitations,  x2  is  in  the  spatial  domain  of  the 

responses  and  t  is  time.  For  linear  operators,  the  Green's  function  G(x,^,t,t)  is  defined 
such  that: 

f-x.i{G(x?,£,t,T)}  =  <S(x,  -  £)5(t  -  t)  (2) 

where  5( )  is  the  Dirac-delta  function.  The  response  of  the  system  can  be  conveniently 
written  in  terms  of  the  Green's  function  as: 

Q(x2,t)  =  JJ  G(x2,£,t,T)w(<j;,T)d£dr  (3) 


since, 
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(4) 


*-,«{Q(x2.t)}  =  JJ  L,,{G(x!,{,!,r))w(^)d(dT 

=  jj  <5(x,  -  £)<5’(t-  T)w(^,T)dgdT 
=  W(x„t) 

Butkovskyi  (1983)  has  proposed  the  introduction  of  a  linear  distributed  block  -  in 
analogy  to  the  lumped  parameter  block  in  classical  control  theory  -  to  represent  the 
input-output  relationship  between  Q(x2,t)  and  w(x1,t).  Thus  the  schematic  of  Figure 
1(a)  is  equivalent  to  the  relationship  expressed  in  Equation  (3  )  .  For  dynamical 
systems  whose  responses  to  stationary  excitations  are  stationary,  i.e  the  Green's 
function  is  stationary  in  time,  the  analysis  can  be  simplified  by  considering  the  Laplace 
transform  of  the  equations  of  motion.  The  role  of  the  Green's  function  is  now  played  by 
the  Transfer  function,  and  the  relationship  of  the  Laplace  transform  of  the  response 
Q(xp,p)  to  that  of  the  excitation  w(x-|  ,p)  is  given  by: 

Q(x2.p)  =  jG(x2f$,p)w($lP)d$  (5) 


where  p  =  o  +  ico  is  the  Laplace  variable;  a  is  the  exponential  growth  rate,  and  co  is 
the  frequency.  This  relationship  is  also  depicted  schematically  in  Figure  1(b). 
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Q(x2,t) 


w(x,.t) 


G(X2>^,t,T) 


(a)  Green's  Function  Representation 


W(x,,p) 


W(x2,|,p) 


Q(x2,p) 


(b)  Transfer  Function  Representation 


Figure  1.  Linear  Distributed  Block. 


Modelling  of  Active  Controller  Attachments 

The  implementation  of  the  active  controller  design  involves  the  application  of 
excitation  forces  at  some  spatial  coordinate  x  =  bp  which  are  proportional  to 

displacements  and  velocities  measured  at  x  =  aj.  For  example,  the  ith  controller 
excitation  force  could  be  written  as: 
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wc(b,,t)  =  g,Q(a,.t)  +  h,Q(a,,t) 


(6) 


where  g,  and  hj  are  the  displacement  and  velocity  feedback  gains  of  the  ith  controller 
respectively,  i  =  1,2,...L,  L  being  the  total  number  of  controllers.  The  Laplace  transform 
of  Equation  (6)  results  in  a  relationship  which  is  used  to  define  the  transfer  function  of 

the  ith  controller  as: 


W0(p)  = 


wc(b„p) 


Q(a„p) 

=  g,+h,P 


(7) 


The  schematic  which  represents  the  combined  interconnected  system  of  the  baseline 
structure  and  the  L  active  controllers  is  shown  in  Figure2.  The  transfer  function  of  the 
combined  system  shown  in  Figure  2  is  given  by  the  following  integral  equation 
(Butkovskyi,  1983): 


W(x,£,p)  =  JWT(x,rjlp)W(rj,^,p)drj  +  W0(x,^,p) 


where, 


WT(x,^p)  =  JW0(xfn.p)£5(77-b1)Wcl(p)5($-a1)d7j 

i.l 

=  iw0(x,b„p)Wo(p)^-a1) 


(8) 


(9) 
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Figure  2.  Schematic  of  Interconnection  of  Linear  Feedback  Controllers 
to  Distributed  Parameter  Baseline  Structure. 

Substituting  Equation(9)  into  Equation  (8)  and  performing  the  integration,  the  result  is, 

W(x,{,p)-£W#(x,b,.p)Wcl(p)W(a..|.p)  +  W0(x,{.p)  (10) 

<  =  1 

In  order  to  solve  for  the  quantities  W(aj,^,p),  i  =  1 ,2 . L,  both  sides  of  Equation  (10)  are 

successively  multiplied  by  5(x  •  am)  and  integrations  are  performed  over  the  x  domain 
for  m  =  1 ,2 . L  to  get: 

W(am.^.p)  =  X'/Vo(a„.b„p)Wtl(p)W(al,e,p)+W0(am,{.p)  (11) 

1=1 

Equation(ll)  is  a  system  of  L  linear  equations  defining  L  unknown  quantities.  If  an 
(LxL)  matrix  [Q]  is  defined  such  that  its  elements  are, 


1 1 


Q,i  =  W0{a(,b„p)WCI(p) 


(12) 


and  if  an  (Lxl)  vector  {«}  is  defined  such  that  its  elements  are, 


a),  =  W(a„£,p) 


(13) 


then  the  system  of  Equation  (1 1 )  for  m  =  1,2 . L  can  be  written  in  a  compact  form  as: 


M  =  [Q]M  +  {y} 


where  {y}  is  an  (Lxl)  vector  whose  elements  are, 


y)  =  W0(a1,^,p) 


From  (14), 

{o)}  =  [l-Q]‘1{y} 

where  I  is  the  (LxL)  identity  matrix. 


(14) 


(15) 


(16) 


Application  of  Computer  Algebra 


The  general  algebraic  form  of  the  transfer  function  of  the  baseline  distributed  system 
is  taken  to  be  (Chen,  1966;  Stakgold,  1979;  Butkovskyi,  1982,  1983;  Keener,  1988; 
Fabunmi,  1989): 


W0(x,£,p)  = 


yf  i  (Po'.ix)(Pok(%) 

p2  "  Pok  J 


(17) 


where  (p0k(x) is  the  kth 


orthonormal  modal  function  for  the  baseline  structure,  p0k  is  the 
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corresponding  modal  parameter  and  f^2  is  the  weighting  factor  or  generalized  mass. 

Although  the  summation  in  Equation  (17)  includes  an  infinite  number  of  terms,  the 
practical  implementation  of  that  expression  requires  that  only  a  finite  number  of  terms 
be  retained.  The  ability  to  retain  a  given  number  of  modes  in  the  algebraic  derivation 
depends  on  the  power  and  memory  of  the  computer  as  well  as  the  number  of  discrete 
modifications  to  the  baseline  structure.  The  substitution  of  Equation  (17)  into  Equation 
(16)  and  the  subsequent  evaluation  and  simplification  of  the  transfer  function  of  the 
combined  system  as  shown  in  Equation  (10)  is  performed  using  the  following  set  of 
MACSYMA  routines: 

WO(EXX,XXSI,PEE):=BLOCK( 

PURPOSE:”EXPRESSION  FOR  TRANSFER  FUNCTION  FOR  BASELINE 
STRUCTURE  -  I.E.  THE  FUNCTION  W0(X,XSI,P)” 

RAT(SUM(‘PHI(EXX,K)*’PHI(XXSI1K)/(PEEA2-(‘P0[K]A2),K,N1  ,N2))/’BSQ)$ 

GAMMA_VECTOR(XXXSI,ARGP):=BLOCK( 

GAMMA  :ZEROMATRIX(NS,  1 ), 

FOR  J  THRU  NS 

DO  SETELMX(W0(‘A[J],XXSI,ARGP),J,1,GAMMA))$ 

OMEGA_MATRIX(ARGP):=BLOCK( 

CAP_OMEGA:ZEROMATRDC(NS,NS), 

FOR  I  THRU  NS 
DO  ( 

FOR  J  THRU  NS 

DO  (W0IJ:W0(‘A[J],’B[I],ARGP), 

SETELMX(WOIJ*WCP[I],I,J,CAP_OMEGA))))$ 
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W(EXX,XXSI,PEE,N11,N22,N33):=BL0CK( 

SCALARMATRIXP:FALSE,N  1  :N1 1  ,N2:N22,NS:N33, 

FOR  KK  FROM  N1  THRU  N2 
DO  STARTP(KK), 

FOR  N  THRU  NS 

DO  WCP[N]:RAT(SUBST(PEE,P,WC[N])), 

OMEGA_  MATRIX(PEE),  GAMMA_VECTOR(XXSI,PEE), 
MATRIX:IDENT(NS)-CAP_OM^oA,OMEGA:ZEROMATRIX(NS,l), 
INVERSE_MATRIX:RAT(ADJOINT(MATRIX))/RAT(DETERMINANT(MATRIX)), 
POLY  :DENOM(INVERSE_MATRIX[  1,1]), 

OMEGA:RAT(INVERSE_MATRIX  .  GAMMA), 
W1:RAT(SUM(WO(EXX,B[I],PEE)*WCP[I]*OMEGA[I,1],I,1,NS)+ 
WO(EXX,XXSI,PEE)),  “DONE”)$ 

W1  gives  the  expression  for  the  transfer  function  of  the  combined  system;  POLY  is  the 
characteristic  polynomial  of  the  combined  system.  In  order  to  cast  the  resultant  transfer 
function  into  the  form  of  Equation(17)  for  the  combined  system,  the  system  parameters 
Plk  are  determined  as  the  roots  of  the  characteristic  polynomial  of  the  system.  In  the 
above  routines,  it  is  possible  to  consider  any  range  of  terms  [Nl  l,N22]  in  the  baseline 
transfer  function  series,  as  well  as  any  number  [N33]  of  discrete  attachments  to  the 
baseline  system.  The  computer-algebraic  results  that  will  be  presented  in  the  next 
section  have  been  generalized  to  the  case  of  an  arbitrary  number  of  discrete 
attachments.  This  is  done  by  mathematical  induction,  based  on  the  results  provided 
by  MACSYMA  for  different  number  of  attachments  specified  in  the  function  calls. 
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V.  COMPUTER-ALGEBRAIC  RESULTS 


Some  results  of  the  derivation  of  the  characteristic  polynomials  for  an  arbitrary  number 
of  attachments  using  one-  and  two  terms  in  the  base-line  transfer  function  series,  are 
presented  in  this  section.  These  derivations  were  performed  on  the  Symbolics  3620.  A 

uniform  one-dimensional  baseline  structure  was  assumed  in  this  study,  so  that  P^2=P2 
for  all  the  k’s. 

One-term  derivation: 

POLYNOMIALS  )  =  C ' 1 '  p2  +  C  (0’ !  (p)  (1 8) 

where, 

C','l(P)  =  -/5Jpk-t(W.(P)«’on,(a,)%n,(b,))  (19) 

1=1 

and, 

cy  >=/?2  (20) 

Note  that  the  polynomial  in  Equation  (18)  is  not  yet  fully  explicit  in  p  until  the 
expressions  for  Wcj(p)  are  substituted  into  Equation  (19).  If  these  functions  are  as 
shown  in  Equation(7)  then, 

c  (p)  =  — /J  JP  On,—  X((9,+  P,P)«>o„,  (a,)«V,(P.)l  (21) 

i=1 

So  that  for  this  class  of  active  controller  design,  and  using  the  one-term  approximation 
to  the  baseline  transfer  function,  the  characteristic  polynomial  in  p  whose  roots  are  the 
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system  parameters  is  given  by: 


POLYNOMIAL^ )  =  {/?2}p2- 


(22) 


Two-term  derivation: 


POLYNOMIAL^)  =  C'V +C  <22!(p)  p2  +  C‘21(P) 


(23) 


where, 


(24) 


C<22,(P )  =  ~  (£  2)  2  £{pOn,  }  "  £  |Wo(P)  X  [  P  Von,  (»,)  <P0n,  (*>,)  ]  J 

=  -(/n2tW-t  jwc(P)L^on,(a,)  ^on,(a,)  J  {j£ 


/?^0n,(b,)' 

(b,)| 


(25) 


and, 
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C'o2’(p)  =  ()32Pon1)(i32pOnJ) 

+x|wci(p) [/!*>„, (aj  />%„,{ a.)  J 


o 

jp<(>  On,(b,)] 

O 

_ i 

Pon,_ 

[/?<?>  0n2(b,)j  / 

L  L 

*11 

i=1  j>i 


W0(p)W,(p)|  |_  <p0„,(a.)  -  «>„,  (a,)  j  | 


|  I  f^On.Cb,)^ 

L^(b.)  "^(b*)  JL0n  (b;) 


(26) 


As  mentioned  earlier,  the  explicit  dependence  of  the  polynomial  on  the  Laplace 
variable  p  will  be  determined  when  the  appropriate  expressions  are  substituted  for  the 
functions  Wcj(p).  If  the  controller  transfer  functions  given  in  Equation(7)  are  substituted 
into  Equation  (25)  and  (26),  then  following  a  collection  of  the  coefficients  of  the  powers 
of  p  in  the  polynomial  of  Equation  (23),  the  result  is: 

POLYNOMIAL^)  =  atp*  +  a3  p3+a2  p2  +  a,  p  +  a0  (27) 

where, 

(28) 
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^TlR,}  I-{aL/JP»>.(a')  J^^b,)}} 

[<P0„,M  -<f0„t(a,)  J  {^,^'jg')}^ 

x  (K,<b,)  -«wwj{£g}) 


l  L 

IS 

1=1  ]>i 


«,  =  1 


h,  [P<P^,M  /*«V,(a.)  J 


Von,  0 


0  pUW.„,<b,)| 


L  L 

-II 

i=1  j>i 


(fl,hi+Olh,)x^  Po n,(a.)  -P«,(a,)  J  {^'(a')}) 
L  «V,<b.)  -?b.,(b,)  J  {**' 


P*.,(b,)l 


and, 


«0  =  (^?Pon,)-(/^2PonJ 

[g,  Ly3<P0n,(a.)  p(p0n,( a,)  J 


+t< 

i  =  1 


Pon 


rp<Po»Sb,Y 


L  L 

-II 

1  =  1  J>1 


L  0  PonJl^^On^b.) 

(g,8,)x(L«V,W  -^,(a.)j{^J|}) 

L  <p°"-(b')  J  { Vcw*(b])} 


(30) 


(31) 


(32) 


For  the  class  of  active  controllers  considered  in  this  study,  i.e  those  with  displacement 
and  velocity  feedback,  the  preceding  expressions  permit  the  explicit  understanding  of 
how  the  feedback  gains  affect  the  resulting  system  parameters.  As  a  matter  of  fact, 
since  polynomials  up  to  the  fourth  order  can  be  solved  in  closed  form  by  computer 
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algebra,  i.  ‘ollows  that  closed  form  expressions  can  be  obtained  for  the  combined 
system  parameters  in  terms  of  the  feedback  gains,  as  well  as  the  location  of  the 
measurement  and  actuation  points.  In  the  following  section,  the  results  for  the  one- 
term  derivation  will  be  used  to  illustrate  a  possible  method  for  determining  the 
controller  feedback  gains  when  it  is  desired  to  achieve  specified  closed-loop  svstem 
parameters.  As  will  be  seen,  this  method  involves  relatively  simple  calculations  which 
can  easily  be  programmed  into  embedded  processors  of  “smart"  structures. 

VI.  CALCULATION  OF  CONTROLLER  FEEDBACK  GAINS 

Consider  a  baseline  structure  with  negligible  damping,  i.e.  p0n1  =  ico0ni  By  virtue  of 
the  weak  coupling  of  the  modes  of  an  undamped  structure,  the  one-term  derivation  is 
adequate  for  estimating  the  closed-loop  system  parameters.  Equating  the  polynomial 
of  Eq.22  to  zero  and  solving  for  the  real  and  imaginary  parts  of  p1n1,  ,  the  frequencies 
and  exponential  growth  rates  of  the  closed-loop  modes  are; 


The  damping  of  the  closed-loop  system  is  controlled  by  the  velocity  feed-back  gains, 
whereas  the  frequency  is  affected  by  both  the  velocity  and  displacement  feedback 
gains.  If  L  velocity  feedback  gains  are  to  be  calculated  directly  based  on  specified 
values  of  the  growth  rates  of  L  system  modes,  Equation  (34)  provides  a  set  of  L 
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equations  for  the  L  desired  values  of  hj.  It  may  also  be  desired  that  there  be  no  shift  in 
the  frequencies  of  these  or  some  other  L  modes.  In  that  case,  a  set  of  L  equations  for 
the  L  values  of  the  displacement  feedback  gains,  gj  ,  can  be  set  up  as  follows: 


£{h,  <?<>,,,  (a.Knjb,)} 


i  =  1 


(35) 


Having  obtained  the  values  of  the  velocity  and  displacement  feedback  gains  in  this 
manner,  it  is  necessary  to  check  the  exponential  growth  rates  of  the  modes  that  were 
not  included  in  the  analysis,  using  Equation  (34).  The  purpose  of  this  check  is  to  verify 
that  there  are  no  modes  for  which  the  exponential  growth  rate  is  positive  -  an 
indication  that  instability  of  that  mode  can  be  induced  by  the  controller. 


Example 

As  an  illustration  of  this  approach,  consider  the  cantilevered  uniform  beam  shown  in 
Figure  3.  The  beam  is  100  meters  long,  with  the  following  cross  sectional  properties: 

flexural  rigidity,  El  =  1.0e8  N-m2  ;  mass  per  unit  length,  pA  =  1  kg/m.  The  natural 

frequencies  and  the  orthonormal  modes  of  the  first  ten  modes  of  this  beam  are  shown 
in  Table  1.  The  active  control  system  consists  of  two  sensors  and  two  actuators.  The 
sensors  are  located  at  spanwise  coordinates  30  and  100  m.  The  actuators  are  located 
at  spanwise  coordinates  40  and  100m.  Consistent  with  the  notation  in  this  paper,  the 
active  control  system  is  made  up  of  four  (4)  controllers:  controller  #1  generates  a  force 
signal  at  x=40  based  on  the  sensor  signal  at  x=30;  controller  #2  generates  a  force 
signal  at  x=l00  based  on  the  sensor  signal  at  x=30;  controller  #3  generates  a  force 
signal  at  x=40  based  on  the  sensor  signal  at  x=100;  and  controller  #4  generates  a 
force  signal  at  x=100  based  on  the  sensor  signal  at  x=100. 


20 


21 


Figure  3.  Example  of  cantilevered  uniform  beam  with  two  actuators  and  two  sensors. 


Table  1.  The  First  Ten  Orthonormal  Modes  of  the  example  Uniform  Cantilevered  Euler-Bernoulli  Beam 

Flexural  Rigidity  EI=1  0e8  N-m2;  Mass  per  unit  length.  pA=l  kg/m;  Length.  L  =  100  m 
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890.732  0.10085  .  0.29845 


Because  of  the  linearity  of  both  the  structural  model  and  the  controller  design,  there 
need  only  be  one  physical  actuator  at  x=40  and  another  one  at  x=100.  The  force 
signals  for  actuators  1  and  3  are  summed  and  applied  to  the  physical  actuator  at  x=40; 
similarly,  the  force  signals  for  actuators  2  and  4  are  summed  up  and  applied  to  the 
physical  actuator  at  x=100.  Also,  only  one  physical  sensor  at  x=30  is  needed  for 
controllers  1  and  2,  and  one  physical  sensor  at  x=100  is  needed  for  controllers  3  and 
4.  The  table  in  Figure  3  specifies  the  sensor  and  actuator  positions  for  the  four 
controllers  in  this  example.  The  objective  of  this  example  is  to  show  how  to  calculate 
the  velocity  and  displacement  feed-back  gains  hj,  gj,  i=1 ,2,3,4  such  that:  (1 )  the  modal 
damping  of  four  selected  modes  are  as  specified  in  advance;  (2)  the  natural 
frequencies  of  four  selected  modes  (not  necessarily  the  same  ones  as  in  (1))  are 
unchanged;  and  (3)  none  of  the  first  ten  modes  is  destabilized  by  the  active  control 
system.  The  restriction  to  the  specification  of  the  parameters  of  oniy  four  modes  is  due 
to  the  fact  that  there  are  only  four  independent  controllers  under  consideration.  The 
feedback  gains  calculated  in  this  manner  are  not  optimal  in  the  usual  sense  of 
minimizing  some  objective  function  which  is  related  to  both  the  response  and  the 
control  power.  The  merit  of  this  approach  lies  in  the  ability  to  concentrate  available 
control  power  on  the  damping  of  certain  modes  which  are  considered  most  responsive 
to  a  given  external  excitation,  without  destablizing  the  other  modes.  Thus  if  the  nature 
of  external  excitation  were  to  change,  and  hence  require  that  some  other  modes  be 
critically  damped,  this  method  affords  a  means  of  adjusting  the  feed-back  gains 
appropriately. 

For  each  of  the  four  modes  for  which  a  desired  damping  ratio  is  specified  in  advance, 
Equation  34  gives  four  linear  equations  for  the  four  unknown  coefficients  hj,  i=1 ,2,3,4; 
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which  are  the  velocity  feed-back  gains.  After  the  velocity  feed-back  gains  have  been 
determined,  Equation  35  is  then  used  to  obtain  the  additional  relations  needed  to 
determine  the  displacement  feed-back  gains  gj ,  i = 1 ,2,3,4;  based  on  the  conservation 
of  the  natural  frequencies  of  the  desired  modes.  For  this  example,  it  is  required  that  the 
first  mode  be  critically  damped,  and  that  the  next  three  modes  be  moderately  damped. 
It  is  also  required  that  there  be  no  shift  in  the  frequencies  of  the  first  four  modes  of  the 
beam.  Table  2(a)  shows  the  specified  damping  ratios  for  the  first  four  modes  of  the 
cantilevered  beam.  For  the  first  mode  to  be  critically  damped,  the  damping  ratio  is 
specified  to  be  unity.  The  results  of  the  calculations  of  the  velocity  and  displacement 
feed-back  gains  are  given  in  Table  2(b)  for  this  example.  Having  obtained  the  values 
of  the  four  pairs  of  displacement  and  velocity  feedback  gains  for  the  four  controllers,  it 
is  now  possible  to  use  Equations  33  and  34  to  recalculate  the  closed  loop  frequencies 
and  exponential  growth  rates  of  all  the  other  modes.  The  results  of  this  calculation  for 
the  first  ten  modes  of  the  example  uniform  cantilever  beam,  are  presented  in  Table 
2(c).  This  example  shows  that  the  feedback  gains  calculated  will  yield  the  desired 
damping  ratios  and  frequency  shifts  without  destabilizing  any  of  the  first  ten  modes  of 
the  example  beam.  So  far,  there  is  no  guarantee  that  all  the  higher  order  modes  will 
be  stable.  If  upon  checking  further,  some  higher  order  mode  is  found  with  a  negative 
damping  ratio,  it  is  prudent  to  perform  the  calculation  again,  including  the  unstable 
mode  among  the  modes  for  which  the  damping  ratio  is  specified.  Further  research  is 
needed  to  develop  a  more  systematic  approach  for  ensuring  the  stability  of  the  modes 
for  which  specific  damping  values  were  not  specified  in  advance. 
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VII.  SENSITIVITY  ANALYSIS 


The  expressions  given  in  Equations  33  and  34  show  that  the  sensor/actuator 
placement  controls  how  the  feedback  gains  influence  the  frequencies  and  damping  of 
the  participating  modes  of  the  baseline  structure.  This  effect  is  more  clearly  displayed 
by  the  sensitivity  of  the  growth  rate  to  the  velocity  feedback  gains  as  well  as  the 
sensitivity  of  the  frequencies  to  the  displacement  feedback  gains: 


<fom, 

<9h, 


2  Bq, 


i 

2(32 


<Pon1(a,)<Pon,(b)) 


(35) 


When  the  sensor  and  actuator  of  the  i-th  controller  are  collocated  (i.e.  a,  =  bj).  the  right 
hand  side  (RHS)  of  Eq.  35  is  positive,  and  all  structural  modes  with  non-zero  values  of 
the  orthornormal  modes  at  the  controller  location,  are  guaranteed  to  be  damped  by  a 
negative  velocity  feedback  gain.  For  controllers  with  non-coltocated  sensor/actuator 
pairs,  the  RHS  of  Eq.35  could  be  positive  or  negative  for  different  modes  of  the 

baseline  structure.  The  baseline  modes  for  which  the  product  <p(aj)cp(bj)  is  positive,  will 
be  damped  for  negative  velocity  feedback  gains,  whereas  those  modes  with  a 
negative  product  of  <p(aj)<p(bj)  could  experience  destabilization  as  a  result  of  negative 
velocity  feedback  gain. 
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<?g, 


(36) 


-  =  0 


=  -j^rZ{h^on,  (a  J<Pon,  (bk )}  j<Pon,  (a>0n,  (b.)  (37) 

As  expected,  the  growth  rate  is  not  affected  by  the  displacement  feedback.  However, 
the  closed  loop  frequency  is  sensitive  to  the  velocity  feedback  as  well  as  the 
displacement  feedback.  For  a  single  sensor/actuator  pair  (L=1 ),  the  sign  of  the  RHS  of 
Eq.37  is  dictated  by  the  sign  of  the  velocity  feedback.  In  the  following  section,  the 
numerical  example  of  a  uniform  cantilevered  Euler-Bernoulli  beam  is  used  to  illustrate 
the  feedback  gain  sensitivities  of  the  closed-loop  growth  rate,  as  the  sensor/actuator 
pair  placement  is  varied,  both  for  collocated  and  non-collocated  controllers. 

Example 

Consider  a  uniform  cantilevered  Euler-Bernoulli  beam,  with  parameters  as  shown  in 
Fig.  4,  where  the  first  five  orthonormal  modes  have  also  been  displayed.  Suppose  a 
single  sensor/actuator  pair  is  to  be  employed  to  control  the  vibrations  of  this  beam 
using  displacement  and  velocity  feedback.  Three  situations  will  be  examined: 

(1)  the  sensor  and  actuator  are  collocated  at  the  same  spanwise  coordinate;  (2)  the 
sensor  and  actuator  are  not  collocated,  with  the  sensor  always  at  the  tip  of  the  beam; 
and  (3)  the  sensor  and  actuator  are  not  collocated,  with  the  sensor  at  midspan.  In 
each  of  these  situations,  the  growth  rate  sensitivity  to  velocity  feedback  gain  for 
different  locations  of  the  actuator,  will  be  studied.  The  growth  rate  sensitivity  to  velocity 
feedback  gain  is  characterized  by  the  partial  derivative  of  the  growth  rate  with  respect 
to  the  velocity  feedback  gain.  As  given  by  Eq.  35,  this  quantity  depends  on  the 
locations  of  the  sensor  and  actuator.  For  vibration  suppression,  it  is  required  that  the 
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Figure  4.  The  First  Five  Orthonormal  Modes  of  a  Uniform  Cantilevered  Euler-Bernoulli  Beam 


growth  rate  be  negative,  and  also  that  its  sensitivity  to  velocity  feedback  gain  be  of  the 
same  sign  for  all  the  modes  of  the  structure.  If  the  sign  of  this  sensitivity  is  not  the 
same  for  all  modes,  velocity  feedback  gains  which  suppress  the  vibration  of  certain 
modes,  will  tend  to  cause  increased  response  others. 

The  plots  of  the  growth  rate  sensitivity  to  velocity  feedback  gains  for  different  locations 
of  the  actuator  are  shown  in  Figures  5-7.  Figure  5  is  for  the  case  where  the  sensor  and 
actuator  pair  are  collocated,  whereas  in  Figures  6  and  7,  the  sensor  and  actuator  are 
not  collocated.  In  Figure  6,  the  sensor  is  always  at  the  tip,  and  Figure  7  is  for  the  case 
where  the  sensor  is  always  at  midspan.  When  the  sensor  is  collocated  with  the 
actuator,  as  in  Figure  5  for  all  locations  of  the  actuator,  or  in  Figure  6  when  the 
actuator  is  at  the  tip,  or  in  Figure  7  when  the  actuator  is  at  midspan,  the  growth  rate 
sensitivity  to  velocity  feedback  retains  the  same  sign  for  all  the  modes  of  the  beam. 
Collocated  sensor/actuator  pairs  near  the  tip  of  the  cantilevered  beam  have  higher 
sensitivities  of  the  growth  rate  of  the  lower  modes  to  velocity  feedback,  than  the  higher 
modes.  This  situation  is  reversed  near  the  root  of  the  beam.  When  the  sensor  is  not 
collocated  with  the  actuator,  there  are  sign  reversals  of  the  growth  rate  sensitivity  to 
velocity  feedback,  for  different  modes  of  the  beam.  In  these  cases,  the  same  velocity 
feedback  gain  which  increases  the  damping  of  certain  modes,  will  reduce  the  damping 
of  others.  These  effects  are  well  known,  and  are  not  surprising.  What  is  interesting  in 
the  present  approach  is  that  these  sensitivities  can  be  quantified  explicitly.  Trade-offs 
can  be  made  with  regard  to  the  selection  of  modes  which  need  maximum  sensitivity,  at 
the  expense  of  other  modes  which  may  not  be  critical  in  a  given  application. 

This  example  has  been  presented  as  a  simple  illustration  of  how  the  results  of  this 
analysis  could  be  used  in  guiding  the  selection  of  sensor/actuator  placements.  In  most 
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Figure  5.  Effects  of  Actuator  Placement  on  Closed-loop  Damping  Sensitivity  to  Velocity  Feedback  -  Collocated  Sensor/Actuator  Pair 
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Figure  7.  Effects  of  Actuator  Placement  on  Closed-loop  Damping  Sensitivity  to  Velocity  Feedback  -  Non  Collocated  Sensor/Actuator 
Pair;  Sensor  @  Midspan. 


practical  applications,  the  orthonormal  modes  may  not  be  available  in  analytical  form. 
It  may  become  necessary  to  conduct  experimental  modal  testing,  in  order  to  obtain  the 
modal  data  required  by  this  approach.  Regardless  of  how  the  modal  data  is  obtained, 
the  simple  expressions  derived  in  this  paper  can  be  used  to  obtain  preliminary 
assessments  of  the  most  advantageous  placement  of  sensor/actuator  pairs  for  active 
vibration  control. 


VIII.  CONCLUSIONS 


This  report  has  presented  results  of  computer-algebraic  derivations  of  the 
characteristic  parameters  of  systems  consisting  of  a  distributed  baseline  structure  and 
output  feedback  linear  active  controllers.  Expressions  which  show  the  explicit 
dependence  of  the  system  parameters  on  the  displacement  and  velocity  feed-back 
gains  as  well  as  the  measurement  and  actuator  coordinates  were  obtained  for  the 
cases  when  the  transfer  function  of  the  baseline  system  has  been  approximated  by 
retaining  one-  and  two-  terms  in  the  infinite  series  which  determine  these  transfer 
functions.  An  immediate  application  of  these  results  is  the  calculation  of  displacement 
and  velocity  feed-back  gains  based  on  requirements  that  certain  closed  loop  modes 
have  specified  damping  ratios.  Because  of  the  simplicity  of  the  calculations  involved  in 
this  process,  it  becomes  practical  to  conceive  embedded  systems  which  permit  the 
"smart”  structure  to  readjust  its  feedback  gains  in  order  to  increase  the  damping  of 
those  of  its  modes  which  are  most  strongly  excited  by  the  external  dynamical  forces. 
The  example  of  a  cantilevered  uniform  beam  was  used  to  illustrate  how  to  implement 
this  methodology  for  a  finite  number  of  controllers.  Further  research  is  indicated  to 
develop  a  systematic  way  of  assuring  that  the  active  control  system  does  not  cause  a 
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destablilzation  of  those  modes  for  which  the  damping  ratios  were  not  specified  in 
advance. 

The  sensitivities  of  closed-loop  system  parameters  to  displacement  and  velocity 
feedback  gains  are  important  considerations  for  the  placement  of  sensors  and 
actuators  for  active  suppression  of  structural  vibrations.  Whereas  optimal  controller 
design  methods  using  modern  control  theory  can  yield  the  optimal  feedback  gains, 
they  do  not  provide  insights  into  to  the  placement  of  the  sensors  and  the  actuators  for 
the  best  effect.  Using  computer  algebra,  explicit  expressions  have  been  derived  for  the 
sensitivities  of  the  closed-loop  system  parameters  such  as  resonant  frequency  and 
exponential  growth  rates,  to  velocity  and  displacement  feedback  gains.  These 
expressions  make  it  possible  to  select  locations  for  sensors  and  actuators,  at  which  the 
closed-loop  system  parameters  have  the  desired  sensitivities  to  the  velocity  and 
displacement  feedback  gains.  Numerical  examples  based  on  a  cantilevered  uniform 
beam  were  used  to  illustrate  the  application  of  this  method,  and  to  show  that  it  is 
possible  to  conduct  a  rational  assessment  of  the  most  advantageous  locations  of 
sensors  and  actuators  for  active  control. 
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